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application in SPH method 
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Abstract — In this paper, we propose a new kernel function 
W e used in the Smoothing Particle Hydrodynamic (SPH) 
method. Three kernel functions are employed to solve the 
differential equations. Numerical results show that W £ has the 
highest efficency. 


Index Terms — kernel function, smoothing particle 
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I. Introduction 

The Smoothing Particle Hydrodynamic method is a 
completely without grid and pure Lagrange particle method 
gradually developed in the nearly 30 years. It is also a 
promising numerical calculations!!] and catch widespread 
concern[2,3] in recent years. 

In the SPH method, the value of function/(x) at the 
point v can be represent by Integration as 

f(x) = lf(x')S(x-x’)dx\ 

Jil 

( 1 ) 

where x e R n ,n = 1,2,3 is the space dimension, £(x) is 
normally Dirac function. 

Replaced S(x — x') by the smooth 

function W (x — x h) , then (1) can be written as 
/(x)*f /(x ')W (x - x h)dx \ 

J Q 

hereW(x — x',/z)is also known as kernel or smooth kernel 

function, h is called the smoothing length. 

According to the writing habits of SPH, the above kernel 
approximate operators can be denoted as 

< f(x) >= J f(x ')W(x - Xh)cbc 

The derivatives of / (x) function is obtained easily by 

using the divergence theorem. It’s worth noting that because 
of the compactly supported of kernel, all derivatives 
of/(x) function is determined by each order derivative of 

the kernel function and the value of / (x). 

In the SPH method, function / (x) is discrete as 
N m 

<f{x)>=Y J — f{x ] )-W(x-x ] ,h). 

mPj 


Where the physical meaning of m. and p ] (j = 1,2,* • *, N) is 

the quality and density of the i-th particle, respectively. 

The kernel of SPH method typically required to meet the 
following conditions: 

1. Regularization: f W(x — x \ h)dx ' = 1; 

J Q 

2. Limit is Dirac function: lim W (x - x h) = 8 (x - x'); 

/*->o v 7 

3. Compact support: W(x — x\h) = 0 when |x - x'| > kh , 
where k is the constant related to the point x; 

4. The distance |x - x'| between x and x' is monotonically 
decreasing and so on. 

Many kernel functions satisfied the above conditions 
have been proposed in the development of SPH method. Such 
as the commonly used Gaussian kernel W G (R , h) (abbreviated 
as W G ), Spline kernel, Index kernel and B-Spline 
kernel W B (R,h) (abbreviated as W B ) and so on. We propose 
a new infinity differentiable kernel function in this paper, and 
compared with W G , W B . To this end, more detailed 

description of these two kernels are introduced in the 
following, for other kernel function can refer[l]. 


The Gaussian kernel W G ( R , h) is defined as[l] 

W G =a d exp(-R 2 ). 

Where the value of a d is related to the space dimension n, i.e. 
n- 1; 


7T m h 

1 


7rh 


TC h 


, n = 2; 


, n- 3. 


Where R = - 


x -x 


Gaussian kernel W G is a better choice when it comes to 
higher derivative because of its Sufficiently smooth, stable 
and high precision computing. However, due to the support 
domain radius of Gaussian kernel functions relatively large, 
therefore, we need a relatively large amount of 
computation 1] for the numerical solution of partial 
differential equations. 
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W B =a d x< 


(2 -Rf 



0<R<1; 

1 < R < 2; 
2<R. 


Where the definition of R is the same as W G , the value of the 
parameters in one-dimensional, two-dimensional and 

1 15 3 

three-dimensional space is — ,-- , - 7 , respectively. 

h 7 nh 1 2jrh 3 

In figure(l), we give function Wc and W B and their respective 
derivative images of the first three in symmetric Interval. 
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(a) W G and its derivative 
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(b) W B and its derivative 
Figure 1 

II. A NEW KERNEL FUNCTION 

Based on level set methods[4]polished to Heaviside 
function, i.e.the smoothness of Heaviside function into 

rr 1 X 1 

H = — arctan — + —. 
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Calculate its derivative 8 S (x) =---- 

n £+ x 

new kernel function is defined in paper[5] 


. Based on this, a 


f 4 1 


W(x-x\s) = < 


n 


l + R 2 ’ 

0, 


1 < R < 1; 

{tf <0}u{tf>l}. 


( 2 ) 

Ijc — jc’l 

Where R =-, h is the parameter, its value determines 

£ 


the smoothness of H e , which affecting the size of the support 


region ofW(x — x\s). function W £ and its derivative images 
of the first three is figure(2). In the following numerical 
examples we found that the value of h for SPH method is also 
affected in solving differential equations accuracy. How to get 
the best value is the content of our further discussion. 
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Figure 2 W e and its derivative 
From the above definition of W(x — x \ s) , it is known 
that this function maintains relatively good properties of 
meshless method. 








wm m » 

— = m s m 




T T 

f \ 

—L 














\ 

1 .. 








T 

1 





/ 




77 

/ \ 




























III. The comparison of three different kernel 

FUNCTIONS FOR SOLVING DIFFERENTIAL EQUATIONS 


In order to further study the properties of the kernel 
function W, We compare SPH method at three kernels, 
i.e. W G , W B , W £ , for the numerical results when solving 
differential equations. 

Given by boundary value problem of second order ordinary 
differential equations: 

d 2 f(x) . 

--— = sin x, 0 <x<n\ 

' dx 2 

f (0) = 0 , /(*-) = 0. 


(3) 

It’s easy to know the analytical solution isf(x) = — sinx. 
According to the SPH method, equation (3)can be 

n ffi 

"(x - Xj , K) = sin x 

J 1 Pi 


Where Xj(j = 1,2,-- -,n) denotes the interval split node, 
parameters K denotes h or £ . 
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n 



w B 

£ 

error 

error 

error 

250 

0.249ft 3 

1.249e-08 

0.0100 

0.0100 

500 

0.124 h 3 

8.869e-09 

0.0050 

0.0050 

1000 

0.062ft 3 

6.285e-09 

0.0025 

0.0025 

2000 

0.031A 3 

4.452e-09 

0.0012 

0.0012 

4000 

0.016ft 3 

3.163e-09 

6.192e-04 

6.192e-04 

8000 

0.008ft 3 

2.284e-09 

3.096e-04 

3.096e-04 


Table 1 SPH method at three kernels, i.e. W G , W B , W £ , the 
compare between error and n of different split interval for the 
solution of boundary value problem of second order ordinary 
differential equations. 

We first consider the impact of parameter £ on the 
numerical solution. In figure(3), we compared the analytical 
solutions and numerical solution of 
parameters £ = 0.249 x h 3 , 

£ — 0.4 x h 3 and £ = 0.2 x h 3 when n = 250 . From the graph 
we can know that numerical solution is very sensitive to 
parameters . How to select the optimal parameters are our 
future research. Table (1) shows the numerical results for 
solving boundary value problems when under three kernels. 
From the table we can see that in the SPH method, by 
selecting appropriate parameters £ under the same split, the 
numerical solution of using weight function W £ has more 
precision than 
W G and W B . 



Figure 4 Take the same subdivision interval, point by point 
error of the SPH method at three kernels 

IV. CONCLUSION 

This paper compares the SPH method for solving boundary 
value problems in three different kernel functions. Numerical 
results show that W £ can get accurate results when select the 
appropriate parameter. Kernel W £ is sensitive to the parameter 
when applied in SPH method, so how to select the optimal 
parameter is our work to be carried out next. 
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Figure 3 SPH method gets numerical solution when take 
different parameters in the kernel W £ 

Figure 4 shows point by point error of the SPH method 
for solving boundary value problem at three kernels. 
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